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Any communication scheme's principal goal is providing error-free data 
transmission. By increasing the rate at which data could be transmitted 
through a channel and maintaining a given error rate, this coding is 
advantageous. The message bits to be transmitted will gradually receive 
more bits thanks to the convolution (channel) encoder. At the receiver end of 
the channel, a Viterbi decoder is utilized in order to extract original message 
sequence from the received data. Widely utilized error correction approaches 
in communication systems for the enhancement of bit error rate (BER) 
performance are Viterbi decoding and convolutional encoding. The Viterbi 
decoder and convolution encoder rate for constraints with lengths of 2 and 6 
and bit rates of 1/2 and 1 are shown in this study in the presence of (1/f) 
noise. The performance regarding the convolutional encoding/hard decision 
Viterbi decoding forward error correction (FEC) method affects the 
simulation outcomes. The findings demonstrate that the BER as function of 


signal to noise ratio (SNR) acquired for uncoded binary phase shift keying 
(BPSK) with the existence of additive white Gaussian noise (AWGN) is 
inferior to that acquired with the use of a hard decision Viterbi decoder. 
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1. INTRODUCTION 

Data transmission through noisy channels is addressed by coding theory, which employs a variety of 
source and channel coding/decoding algorithms. Forward error correction (FEC) schemes frequently employ 
convolutional coding. The widely utilized Viterbi algorithm serves as the foundation for the primary 
decoding method for convolution codes. Convolutional encoding with the Viterbi decoding represents potent 
FEC approach which is especially well-suited to a channel where transmitted signal is mostly contaminated 
by additive white Gaussian noise (AWGN) in conventional systems [1]. However, in fact, colored noise is 
present [2]. Through carefully designing redundant information to be added to data that is being transmitted 
via the channel, FEC aims to increase the channel's capacity [3]. For the most part, Viterbi algorithm 
conducts maximal likelihood decoding to fix errors in received data that is brought on by channel noise. 
Therefore, reduce bit error rate (BER) to enhance performance. Because Viterbi decoding has set decoding 
time, it’s advantageous for the implementation of hardware decoder. The specifications for Viterbi decoder, a 
processor that executes the Viterbi algorithm, vary depending on the application. Viterbi algorithm uses the 
greatest resources and is the most effective [4]-[6]. In Figure 1, a general block architecture regarding a 
digital communication system with the use of convolution encoder and Viterbi decoder has been shown in the 
presence of colored noise. 
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Figure 1. Diagram of the digital communication system utilizing convolution encoder and Viterbi decoder in 
the presence of the colored noise 


2. METHOD 

Shift registers are used in fixed numbers to create convolutional encoders. The encoder output is 
created through adding bits in shift register, which each input bit enters. The quantity of modulo 2-adders that 
are utilized in combination with shift registers specifies the number of the output bits [1], [4], [7]. 


2.1. Encoder parameters 

Three parameters are frequently used to describe convolutional codes: n, k, and m, in which n 
represent number of the output bits, m represent number of the memory registers and k represent number of 
the input bits. The code rate, which has been denoted by k/n, is a measurement of a code's bandwidth 
effectiveness. Except for the applications of deep space, when code rate values as low as 1/100 or longer 
could be used, k and n parameters typically vary between 1 and 8, code rate between 1/8 and 7/8, and m 
between 2 and 10. The parameters (n, k, and L) could also be used to specify convolutional codes. L is 
referred to as constraint length of code and represents the amount of bits in an encoder memory which have 
an impact on how the n output bits are generated. Instead of (n, k, and m) codes, convolution codes that have 
been presented here will be denoted as (n, k, and L) codes [6], [7]. 


2.2. Structure of the encoder 

The connections of the encoder have been described by term generator polynomial (g). The 
generator polynomial for that output bit refers to decision of which bits (in memory registers) should be 
added (with the use of modulo-q adders) for producing the output bits. For any m order code, there are a 
variety of polynomials to choose from. In most cases, using computer simulations, good polynomials are 
discovered by trial and error. The next presumptions were established for better understanding a 
convolutional encoder's working and the FEC method [8], [9]: i) A (2, 1, 3) convolutional encoder has been 
utilized, ii) a three-bit sequence of input has been utilized by bits [1 0 1], and iii) 2 generator polynomials 
have been utilized, which were specified by bits [1 0 1] and [1 1 1]. 

Constructing convolutional encoder is simple. Initially, m boxes are drawn in order to represent m 
memory registers. Next, the n output bits by n modulo-2 adders are drawn. Lastly, we use the bits defining 
generator polynomials for connecting the memory registers to adders. A (2, 1, 3) convolutional encoder has 
been depicted in Figure 2. The 2 generator polynomials that have been designated by bits [1 O 1] and [1 1 1] 
will be utilized for encoding the three-bit input sequence [1 0 1] using this encoder. The output bits 1 and 2 
are represented by v1 and v2, respectively, whereas ul represents the input bit. The memory registers’ initial 
states, initially set to the value of 0, are represented by u-1 and u-0. 


2.3. Convolutional representation 

The convolutional encoder's logical operation was shown in the preceding section. A look-up table, 
also known as table of state transitions, could be used by the convolutional encoder to do the encoding. There 
are four items in the state transition table: i) input bit, ii) encoder state that is 1 of 4 potential states 
(00 01 10 11) for (2, 1, 3) convolutional encoder, iii) output bits that for (2, 1, 3) convolution encoder are: 
00 01 10 11, due to the fact that only 2 bits would be output, and iv) output state which represents the 
following bit’s input state. The output state which is state transition table for (2, 1, 3) convolutional encoder 
is as it has been clarified in Table 1. There are graphically 3 ways for the representation of encoder for the 
purpose of gaining a better understanding concerning its operation [10], [11]. Those include: i) tree diagram 
representation and ii) state diagram representation. 
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Figure 2. A (2, 1, 3) convolutional encoder 


Table 1. State transitions for (2, 1, 3) convolutional encoder 


Input bit Input state Output bits Output state 
0 00 00 00 
1 00 11 10 
0 01 11 00 
1 01 00 10 
0 10 10 01 
1 10 01 11 
0 11 01 01 
1 11 10 11 


2.3.1. Representation of the state diagram 

The state diagram for (2, 1, 3) convolution encoder is as has been depicted by Figure 3. Contents of 
two shift register cells have been referred to as a state. A state is represented by each circle. The encoder can 
be found in one of such states at any one time. As bits arrive, lines from and to it display possible state 
transitions. The arrival of bit | or the arrival of bit 0 are the only 2 cases that can occur simultaneously. The 
encoder is able to "jump" into a different state as a result of each of these two events. As can be seen, the 
state diagram contains the same data as the state transition table but does it graphically. The dotted lines 
exhibit arrival of bit 1, while solid ones show arrival of bit 0. Each case's output bits are displayed on a line, 
and an arrow points to the state transition [12], [13]. 


Figure 3. Diagram for (2, 1, 3) convolutional encoder 


2.3.2. Tree diagram representation 

Figure 4 displays a tree diagram for the convolutional encoder (2, 1, 3). State diagram entirely 
describes the encoder, however it’s difficult to utilize it to monitor encoder transitions over time. Due to the 
fact that the diagram does not depict time history, this is obvious. The tree diagram gives the state diagram a 
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time dimension by depicting time passage as we delve farther to the 3 branches. Which makes it superior 
method for describing convolutional codes compared to the state diagram. With the tree diagram, branches of 
a tree are travelled based upon arrival of bit 1 or bit 0, rather than "jumping" from one state to another. In the 
encoder's input, if bit 0 is received, we move up the tree, and if bit 1, we move down. The output state is 
represented by the bits in parenthesis after the first two bits, which stand for the output bits [14], [15]. 
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Figure 4. Tree diagram for (2, 1, 3) convolutional encoder 


3. VITERBI DECODING 

The process of the recovery of encoded input data stream at receiver after it was transmitted over a 
channel has been referred to as the channel decoding. Sequential decoding and maximum likelihood 
decoding, which is usually referred to as the Viterbi decoding, represent the 2 main channel decoding types 
for convolutional codes. Maximum likelihood decoding, which is often known as Viterbi decoding, has been 
created by Viterbi, a founding member of Qualcomm corporation. Error bounds for the convolutional codes 
and asymptotically optimal decoding algorithm is his seminal study on the method. Viterbi decoding is more 
suitable for hardware implementation because, unlike sequential decoding, decoding time is fixed rather than 
variable. Here, we methodically eliminate alternatives at each trellis level [16], [17]. The following criteria 
were utilized in order to narrow the options: 
- Errors happen in an infrequent manner. Error likelihood has been found small 
- The likelihood of 2 consecutive errors is considerably smaller when compared to one error that is the 

errors are distributed in a random manner 

The entire received sequence of some specific length is examined by Viterbi decoder, which 
calculates a metric for every one of the paths and bases its decisions on that metric. Up until 2 paths meet at 
one node, all paths are taken. One of the two paths is subsequently chosen depending on an earlier decision. 
The chosen paths are referred to as the survivors. There are a total of 2N potential received sequences for an 
N-bit sequence. Only 2kL of these are valid. Just one survivor path is selected depending on a decision when 
two paths converge on a single node. There are two ways to make this decision, leading to the 2 Viterbi 
decoding forms [18]: i) hard and ii) soft decision Viterbi decoding. 


3.1. Trellis explanation 

The hard decision Viterbi decoder's actual working is depicted in the following figure. Depending 
on the example taken, a trellis is drawn for each time tick of decoded three-bit input stream [1 0 1]. The hard 
decision Viterbi decoder assumes that corrupted data bit stream at input is [01 10 00 10 11]. The path that has 
the maximum metric has been sought out and a winner path is traced after five stages. The decoded sequence 
is represented in Figure 5 and follows the path made up of states 00, 10, 01, 10, 01, 00, which correspond to 
the bit values 10100. It can be seen the way that the hard decision Viterbi decoder generates decoded data bit 
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stream from convolutional encoded input data bit stream transmitted through a coloured noisy channel from 
transmitter by applying maximum Hamming distances [19]. 
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Figure 5. Decoded sequence 10100 for noisy encoded bit stream 01 10 00 10 11 


4. COLORED NOISE GENERATION USING MATLAB 
In various applications, white noise with a Gaussian distribution is assumed. The definitions of the 


autocorrelation and power spectrum are [12]: 


reelm] = E {elmJefm + kl} = 026[m] a) 
See [e21] = F.T. yor {r..[k]}} = Noms <f< S a 


The power spectrum density (PSD) of the white Gaussian noise represents constant value over full 
range of the frequencies, all of the frequency values range with no/2 magnitude. For some certain time 
instant, it had shaped the probability distribution function pdf p,(€) that has been given as (3) [20]: 


_(€-He)? 


e 2o (3) 


pe(€) = dat 


in which the mean value is represented by ue and the standard deviation by oe. Since all of the samples are 
Gaussian and share same statistical characteristics, delta function on autocorrelation functions indicates that 
nearby samples are not dependent. The observed samples have been thought to be independently identically 
distributed (i.i.d.). Figures 6 and 7 display the AWGN's time representation, PSD, respectively [21]. Figure 8 
shows the autocorrelation function for AWGN and its clear similar to delta function. 
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Figure 6. Time representation of the AWGN Figure 7. Power spectral density of the AWGN 
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Figure 8. Autocorrelation function of the AWGN 


According to Figure 9, white Gaussian noise is processed through a low-pass filter to produce color 
noise. It is coloured with 1st order low-pass Butterworth filter with cut-off frequency of 0.1 fs and added to 
signal of interest. This filter represents all-pole filter. The equation: gives its magnitude squared frequency 
response as (4) [22], [23]: 


|H(w)|? = (4) 


ae 
1+(“/y) 


in which w represent the digital frequency, w (c) represent the 3 dB cut-off digital frequency, N represent the 
filter order. Consequently, the frequency response of a first order Butterworth filter is proportional to (5): 


|Hw)P as (5) 


White 
Gaussian Noise 


Figure 9. Generation of colored noise from white noise 


Figure 10 shows the time representation of generated colored noise obviously the difference with the 
Figure 6. The PSD regarding the generated colored noise is shown in Figure 11; it is evident that the noise 
isn’t constant across the entire range of the frequencies and that it is proportional to 1/f noise. The 
autocorrelation function related to the colored noise is illustrated in Figure 12. The autocorrelation functions' 
sinc function indicates that nearby samples are dependant when all of the samples are Gaussian and have the 
same statistical characteristics [24], [25]. 
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Figure 10. Time representation of the colored noise Figure 11. PSD of the colored noise 
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Figure 12. Autocorrelation function of the colored noise 


5. RESULTS 

The performance regarding the convolutional encoding/hard decision Viterbi decoding FEC method 
affects the simulation results. Those performance factors aid in the investigation of how the FEC method 
might be used in practical settings. The next performance factors are taken into account: i) encoder memory 
size and ii) signal to noise ratio (SNR). A hard decision Viterbi decoder with 1/2 code rate and memory m=3 
is used to display the BER as an SNR function in Figure 13. SNR (Eb/No) in dB obtained with our 
MATLAB implementation utilizing the collection of common polynomials is shown in Table 2. 


Plot of BER vs. Eb / No 


Uncoded BPSK AWGN 
—4— code rate = 1/2 , m = 2 , hard-decision, colored noise 


Bit Error Rate (BER) 


Figure 13. BER plot versus Eb/No for a rate 1/2, m=3 encoder, with the consideration of hard decision 
Viterbi decoder 


Table 2. SNR (Eb/No) in dB for the rate 1/2, m=3 encoder 


BER Uncoded BPSK Coded 
107? 3.5 3.5 
103 7.1 6 
10% 8.3 73 


The performance of a hard decision Viterbi decoder with a 1/3 code rate and memory m=3 is 
depicted in Figure 14. The SNR (Eb/No) in dB obtained with the proposed MATLAB implementation 
utilizing common polynomials is shown in Table 3. The performance related to a hard decision Viterbi 
decoder with a 1/2 code rate and memory m=6 is depicted in Figure 15. SNR (Eb/No) in dB that had been 
obtained with our MATLAB implementation utilizing common polynomials is shown in Table 4. 
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Figure 14. Plot of BER versus Eb/No for a rate 1/3, m=3 encoder, taking under consideration a hard decision 
Viterbi decoder 


Table 3. SNR (Eb/No) in dB for rate of 1/3, m=3 encoder 


BER UnCoded BPSK Coded 
10° 7 4.8 
10% 7.4 6 
10° 9.2 6.8 


A Plot of BER vs. Eb / No 
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Figure 15. BER Plot vs. Eb/No for a rate 1/2, m=6 encoder, taking under consideration a hard decision 
Viterbi decoder 


Table 4. SNR (Eb/No) in dB for 1/2 rate, m=6 encoder 


BER UnCoded BPSK Coded 
103 7 4.2 
10% 8.3 6.8 


6. CONCLUSION 


The design of a Viterbi decoder and a convolutional encoder, which could encode a digital data bit 
stream and produce a code word which could be transmitted and decoded at the destination. For each 
constraint length, rates of 1/2 and 1/3 were used when designing the encoder, which had constraints of 
lengths 2 and 6. The Viterbi decoder has been engineered for the determination of decoding path with least 
design of metrics to be sent to the port of the decoder output. Design of a Viterbi decoder and convolutional 
encoder with the use of MATLAB has proven successful, with results in terms of BER vs SNR in the 
presence of colored noise. Future work could apply Turbo codes and Reed-Solomon coding for the signal 
corrupted by colored noise using the FEC approach. 
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